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Abstract. Following the Flory ideality hypothesis intrachain and interchain excluded volume interactions 
are supposed to compensate each other in dense polymer systems. Multi-chain effects should thus be 
neglected and polymer conformations may be understood from simple phantom chain models. Here we 
provide evidence against this phantom chain, mean-field picture. We analyze numerically and theoretically 
the static correlation function of the Rouse modes. Our numerical results are obtained from computer 
simulations of two coarse-grained polymer models for which the strength of the monomer repulsion can 
be varied, from full excluded volume ('hard monomers') to no excluded volume ('phantom chains'). For 
nonvanishing excluded volume we find the simulated correlation function of the Rouse modes to deviate 
markedly from the predictions of phantom chain models. This demonstrates that there are nonnegligible 
correlations along the chains in a melt. These correlations can be taken into account by perturbation theory. 
Our simulation results are in good agreement with these new theoretical predictions. 

PACS. 61.25.H- Macromolecular and polymers solutions; polymer melts - 61.20.Ja Computer simulation 
of liquid structure 



1 Introduction 

In 1953 P.E. Rouse [1] proposed a model to describe the 
dynamics of a polymer chain, which has become an impor- 
tant concept in polymer physics [2j[3]. The Rouse model 
assumes the chain to be a sequence of Brownian beads. 
The beads are connected by entropic springs and immersed 
in a structureless medium which exerts on every bead two 
forces, a local random force and a local friction force. Both 
forces are linked by the fluctuation-dissipation theorem to 
ensure correct equilibrium properties [2]. 

This bead-spring model thus considers only local inter- 
actions. Nonlocal interactions, such as hydrodynamic or 
excluded-volume forces between distant beads along the 
chain, are ignored. While this assumption is certainly not 
appropriate for dilute solution in good solvents, it may be 
valid in concentrated solutions or polymer melts, where 
both nonlocal interactions are expected to be screened [21 
|3J. In particular in polymer melts, the screening is sup- 
posed to extend down to the monomer level, implying that 
polymer conformations correspond to those of ideal ran- 
dom walks. Hence, it is generally believed that the Rouse 
theory provides, at long times and large length scales, a 
viable description of the conformational dynamics of poly- 
mer melts if entanglements with other chains, giving rise 
to reptation motion [SEIH], are not important. 

Experimental or computational scrutiny of the Rouse 
model has thus focused on the behavior of short chains in 
a melt (for review see e.g. [HE]). These tests reveal that 



the Rouse model represents a good approximation, but 
quantitative agreement is hard to obtain. To explain the 
observed deviations there are, roughly speaking, two main 
ideas in the literature: additional intrachain contributions 
not accounted for by the theory (chain stiffness [§l[7l[8], 
local excluded volume effects [MTU], torsional transitions 
[rilfT2l[T3] ) and multi-chain effects which invalidate the 
phantom-chain-in-a-structureless-medium approach of the 
Rouse model [5]. 

In the present article, we provide further evidence for 
the latter point of view. However, there is one important 
respect in which the objective of our discussion is limited. 
The Rouse model is, in the first place, an attempt to de- 
scribe dynamic features of polymer melts. Here we shall 
not be concerned with dynamics, but rather study, by the- 
ory and simulation, the initial (static) value of the corre- 
lation function of the Rouse modes. This initial value in- 
forms us about conformational properties of a chain in the 
melt, and clearly, an understanding of these equilibrium 
features is a prerequisite for an extension of the discussion 
to polymer dynamics. Our key point is that the conforma- 
tional properties, even for very flexible chains, are more 
complex than commonly assumed because chains are not 
phantom chains, free of any interaction. In the melt, poly- 
mer segments of size s (s> 1) experience an effective re- 
pulsion resulting from chain connectivity and incompress- 
ibility of the melt. This repulsion entails systematic de- 
viations from 'ideal' random-walk-like conformations and 
thus leads to a violation of the 'Flory ideality hypothesis' 
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|14j . i.e. of a central concept of modern polymer physics. 
The consequences of these deviations have recently been 
explored in real (intrachain correlations [151116] ) and recip- 
rocal space (form factor [171IT8] , collective structure factor 
(H)]). The purpose of the present work is to extend this 
discussion to the static Rouse modes. 

The outline of the article is as follows. We begin by pre- 
senting the simulation models and techniques used in this 
study (Section [2]). Section [3] compares predictions from 
phantom chain calculations to simulation results of poly- 
mer melts with full excluded volume interaction between 
the monomers. The comparison reveals systematic devia- 
tions between theory and simulation for the static corre- 
lation function of the Rouse modes. These deviations can 
be rationalized by a theory which accounts for corrections 
to chain ideality in polymer melts (Section F4|. Section [5] 
extends the previous discussion to the situation where the 
phantom chain limit is approached by gradually making 
the monomer- monomer interaction softer. Also this case 
can be understood by taking corrections to chain ideality 
into account. The final section (Section[6]) presents a brief 
synopsis of our results. 
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2 Simulation models and techniques 

By computer simulations we examine polymer melts of two 
coarse-grained models |20], a bead-spring model and the 
bond fluctuation model. For both models we study a ver- 
sion which precludes monomer overlap ('hard monomers') 
and another one with finite energy penalty for monomer 
overlap ('soft monomers'). The latter version allows us 
to switch gradually from full monomer excluded volume 
interactions, the standard choice in simulations, to phan- 
tom chain behavior, the situation considered by the Rouse 
model. 

In the following we briefly describe the simulation mod- 
els and techniques. Further details may be found in Refs. 

HUE]. 



2.1 Bead-spring model 



Fig. 1. Illustration of the modeling of soft monomers for the 
BSM [panel (a)] and the BFM [panel (b)]. For the BSM a 
force capping with parameter f c = 1.781 is shown. This gives 
r c = 1.0535 and an energy penalty for monomer overlap of 
E = U Ii b{r c )+f c r Q = 1.94. For the BFM every doubly occupied 
lattice site is penalized by an energy cost of E/8 so that the 
energy penalty for full overlap of two monomers equals E. 



C = 4eo/27. Intrachain and interchain interactions of non- 
bonded monomers are thus purely repulsive. The param- 
eters of the bond potential are adjusted such that the av- 
erage bond length I w l\, = 0.97234 <7o at the monomer 
density p = 0.84 a^ 3 is very close to that of the standard 
Krcmer-Grest model [2"6ll27j . In the following, we report 
all data in reduced units, that is, energies are measured in 
units of k^T (Boltzmann constant /cb = 1) and lengths in 
units of the monomer diameter do- 



2.1.1 Excluded volume chains 

The bead-spring model (BSM) is derived from a model 
employed in simulations of polymer crystallization [22,23, 
24 ,25]. It is characterized by two potentials: a harmonic 
bond potential, 



U h (r) = \k h {r-l h f 



(1) 



where k h = 535.46 k^/T /a^ and l h = 0.97234 a , and a 
nonbonded 9-6 Lennard- Jones (LJ) potential 



U ab (r) = 




+ C r <r- 



mm 3 



(2) 



with Sq = 1.511fceT. The LJ potential is truncated at 
its minimum, r m j n = (3/2) 1 / 3 cro, and shifted to zero by 



2.1.2 From hard to soft monomers 

Equations (JTJ) and ensure that chains cannot cross each 
other and monomers are impenetrable ('hard monomers'). 
We 'soften' these constraints by introducing a force cap- 
ping through the following modification of the nonbonded 
interactions 



M - fU nh (r c ) + f c (r c -r) 
(rJ ~\U nh (r) 



< r < r r 



r > r c 



(3) 



Here the 'capping' distance r c is defined by the condition 
that the nonbonded force at r c has the prescribed value / c , 
i.e. — dU n b(r)/dr\ r=ra = f c . This definition implies that 
r c tends to for / c 3> 1 and to r min for f c — ► 0. (See 
Fig. HJa) for an illustration.) For distances smaller than 
r c Eq. ([3]) thus replaces the steep rise of U n b(r) by a much 
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weaker linear increase, entailing a finite energy penalty 
E = U n b(r c ) + f c r c for two overlapping monomers (r = 0). 
The smaller / c , the softer the monomers. In particular, if 
f c = 0, we recover phantom chain behavior, i.e. U^(r) = 
0, because r c then equals r m i n and U n h{r) vanishes for 

In this way, we also simulate two phantom chain mod- 
els, a (generalized) freely-jointed chain (FJC) model char- 
acterized only by the bond potential of Eq. (UJ), and a 
(generalized) freely-rotating chain (FRC) model which, in 
addition to Eq. JTJ), also presents a potential for the bond 
angle 6, 

U ang (8) = fc e [l-cos2(0-0 o )] , (4) 

with kg = 10 and 0q = 109°. For this potential we find 
that the average cosine of the bond angle is a = (cos 0) = 
-0.320713 w cos O . 



are carried out by a mixture of local, slithering-snake, and 
double-bridging MC moves which allow us to equilibrate 
polymer melts with chain lengths up to N = 8192 [To] . 

2.2.2 From hard to soft monomers 

In analogy to the BSM we also study soft monomers by 
introducing a finite energy penalty E/8 for a doubly occu- 
pied lattice site (see Fig.[JJ). This implies that full overlap 
between two monomers leads to an energy cost of E (as 
for the BSM, cf. Section 12. 1.2|) . A local or slithering-snake 
move, leading to N ov double-occupancies, gives rise to a 
total energy Ti/k-QT = N ov E/8. With the energies of the 
final (7if) and initial configurations (Tii) we accept the 
move according to the Metropolis criterion [2"0ll3"T] with 
probability min(l, exp[— (Tii — TL{) / k^T]) . 



2.1.3 Simulation aspects 

With these models we perform molecular dynamics (MD) 
simulations at constant temperature T = 1 (Langevin 
thermostat with friction constant 7 = 0.5) and constant 
monomer density p = 0.84, the typical melt density of 
the Kremer-Grest model [261127]. The equations of mo- 
tion are integrated by the Velocity- Verlet algorithm [2"5] . 
For 'hard' monomers we combine the MD with double- 
bridging Monte Carlo (MC) moves [201I27] to speed up 
the decorrelation of large-scale conformational features. 
As only few of these MC moves are accepted per unit 
time, this does not deteriorate the stability or accuracy of 
the MD. The MC moves, however, considerably improve 
the statistics for large chain lengths. For 'soft' monomers 
only MD is used. We will discuss data for chain lengths 
(number of monomers per chain) N = 64, 256, 512, and 
1024, obtained from simulations of periodic systems of lin- 
ear size L < 62. For p = 0.84 these systems contain up to 
196 608 monomers. 



2.2 Bond fluctuation model 

2.2.1 Excluded volume chains 

We also examine the three-dimensional bond fluctuation 
model (BFM) on a cubic lattice [29l|30|. Each monomer 
occupies a cube of eight adjacent sites, the length of the 
bonds between connected monomers along a chain are 
allowed to fluctuate in the range from 2 to \/l0 lattice 
constants (the lattice constant will be the length unit in 
the following), and double occupancy of lattice sites is 
forbidden by a hard-core interaction between monomers. 
The system is athcrmal, the only control parameter be- 
ing the monomer density p. Melt conditions are realized 
for p = 0.5/8, where half of the lattice sites are occupied 
[30] . We use periodic simulation boxes of linear dimension 
L = 256 which contain pL 3 « 10 6 monomers. These large 
systems eliminate finite-size effects, even for the longest 
chain lengths studied (256 < N < 8192). The simulations 



3 Rouse model: predictions and observed 
deviations 

The basic variables of the Rouse model are the Rouse 
modes X p . Here we introduce two definitions for X p , de- 
pending on whether we consider a discrete (simulation) 
model or a continuous (theoretical) model. 

For the discrete case let r„ be the position of monomer 
71 (re = 1, ... , N). The Rouse modes are defined by [32] 



1 " (n - l/2)jwr 

71=1 



N 



(5) 



with p = 0, . . . , N — 1. In the continuum limit we will use 
the notation of Ref. [2] and write 

1 f N n ir 

X p = -y o dr l r(n)cos^ (p = 0, 1, 2, . . .) . (6) 

Here the monomer index n is a continuous variable ranging 
from to the total number of bonds N. 

Our analysis focuses on the static correlation functions 
of the Rouse modes p and q. In the continuum limit these 
functions are given by 



pi 



(X p ■ Xg) 

rN r N 



1 / f Tl 7T 

= — k I dn dm (r(n) ■ r(m)) cos cos 

N Jo Jo ' N 



npn mqTT 



2N 2 



N r N 



dn J dm ([r(n) — r(m)] 2 ) x 

npn mqn 



cos cos ■ 

N N 



(7) 



In dense melts it is commonly assumed that intrachain and 
interchain excluded volume interactions compensate each 
other down to the scale of a monomer [5] so that a mean- 
field picture should apply. Static equilibrium features can 
be obtained by treating the polymer liquid as an ensem- 
ble of independent chains displaying ideal random-walk- 
like conformations. If this was true, single-chain models, 
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Fig. 2. Correlation function of the Rouse modes C pp versus 
p/N for the BSM (left ordinate) and the BFM (right ordinate). 
All simulation results shown refer to chains with full excluded 
volume between the monomers. Following Eq. JS} the ordi- 
nate is scaled by 8N/b% (BSM: b c = 1.338, BFM: b c = 3.244 
|16]). The solid lines indicate the predictions of the FJC model 
[Eq. ©]. 



i.e. phantom chains, should suffice to fully describe chain 
conformations. For (some) phantom chain models, the cor- 
relation function of the discrete Rouse modes can be calcu- 
lated. Therefore, it should be possible to predict, say, the 
p-dependence of the diagonal elements of the simulated 
Rouse mode matrix, C pp , from an appropriate phantom 
chain model. In the following we want to provide evidence 
against this expectation. 

The simplest phantom chain model is a freely-jointed 
chain (FJC) model. For the FJC model the result for C pq , 
obtained with the discrete Rouse modes [Eq. J5J], reads 
(see e.g. [35]) 



1 



sin(p7r/2iV) 



p/at«i 



1 Nbl 



2tt 2 



V 



(for p > 0) . (8) 



The limit p/N <C 1 coincides with the result obtained for 
the continuum model [2] . In Eq. ([8]) b c denotes the effective 
bond length defined by b% = R^/N from the end-to-end 
distance R e of asymptotically long chains [2] . 

Figure [2] compares Eq. ([8]) to simulation data for ex- 
cluded volume chains of the BSM and the BFM. The data 
for all chain lengths collapse onto a common curve which 
appears to agree well with [sin(p7r/2iV)]~ 2 , if p/N < 0.01. 
(In fact, the agreement is not as good as it seems; we 
return to this point below.) For larger p/N, however, de- 
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Fig. 3. Correlation function of the Rouse modes C vv ver- 
sus p/N for BSM phantom chains (FJCs and FRCs, cf. Sec- 
tion [2X2j| of lengths N = 256 and TV = 1024. C pp is scaled 
by 2(irp) 2 /N so that the ordinate tends to be for p/N — > 0. 
The simulation results for FJCs and FRCs are compared to 
Eqs. |(8j| and (|9]), respectively. These comparisons utilize the 
following parameters: for FJCs, bl = 0.9532 and for FRCs, 
tit = Coo x 0.9532, where the characteristic ratio is given by 
Coo = (1 - a) /(I + a) (with a = -0.320713). The values of b c 
and a are determined in the simulation. 



viations occur. The FJC model overestimates the corre- 
lation, especially for p/N > 0.1. In [9] and more recently 
also in [TU] it was argued that these modes are dominated 
by the microstructure of the simulated chain model. For 
instance, p/N > 0.3 corresponds to subunits of a trimcr 
and smaller. On these local scales, the fact that the re- 
pulsive monomer interactions of the BSM and the BFM 
avoid immediate backfolding of the chain and thus con- 
fer some intrinsic stiffness to the polymer should be taken 
into account. 

The simplest way to achieve this consists in replacing 
the FJC by a (generalized) freely- rotating chain (FRC). 
One can introduce the bond angle as a further degree of 
freedom and still carry out the summation to determine 



Cpp. For N > 1 the result reads [9] 



8N 

bl pp 



sin(p7r/2iV) 



4a 



2a cos(pn/N) 



(9) 



where a = (cos 6) . 

Figure [3] shows that Eqs. © and © agree respectively 
with the simulation results for the corresponding phantom 
chain models. However, they do not for the full-excluded- 
volumc chains in the melt. This is demonstrated in Fig. [5] 
for the FJC model and in Fig. H for the FRC model. To 
emphasize the differences between the numerical results 
and FRC prediction the simulation data are divided by the 
asymptotic behavior expected (and found) for p/N <C 1 
[Eq. ©]. With increasing p/N, the simulation data show 
a continuous depression, which is not reproduced by the 
FRC model. (The FJC model has not even a minimum; 
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Fig. 4. Upper panel: Correlation function of the Rouse modes 
C pp versus p/N for BSM excluded volume chains. C pp is scaled 
by 2(irp) 2 /Nbc so that the ordinate tends to 1 for p/N — > 
(b c = 1.338 [E]). The solid line presents Eq. (|19p with param- 
eters from Table [1] The dashed line indicates the prediction of 
the FRC model, Eq. JjJ), with a = —0.1948 (measured in the 
simulation). Lower panel: Same as in the upper panel, but for 
BFM excluded volume chains with b c = 3.244 (Table [TJ) and 
a = —0.1055 (measured in the simulation). 



cf. Fig.[5J) This depression is, however, compatible with a 
theory that accounts for the impact of residual excluded 
volume interactions in the melt. We will sketch this theory 
in the next section. 



4 Static Rouse mode correlations: corrections 
to chain ideality 

By means of the random-phase approximation (RPA) Ed- 
wards derived the pair potential between two monomers in 
a dense, three-dimensional multi-chain system of asymp- 
totically long chains [21155], 



v(r) 



5(r) 



cxp(-r/£) 
4-7rr£ 2 



(10) 



This potential consists of two terms. The first term repre- 
sents the bare repulsive interaction between two monomers. 
It is of very short range and characterized by the excluded 



volume parameter v. The latter is related to the compress- 
ibility of the multi-chain system [see Eq. (|24p ] . The second 
term in Eq. (|10[) results from the repulsion of all monomers 
in the system whose compound effect is to attenuate the 
bare interaction between the tagged monomer pair. It is of 
range £, where £ = b/(12pv) 1 / 2 is the (Edwards) correla- 
tion length of collective density fluctuations, and b denotes 
the effective bond length of an ideal polymer chain with 
all interactions switched off (v = 0). 

In a dense (three-dimensional [34]) system we expect 
both terms to nearly compensate each other and v(r) to 
be small. This suggests that a first-order perturbation cal- 
culation should be appropriate to explore the influence of 
v(r) on the conformational properties of a polymer melt. 
Following Eq. Q C pq can be obtained from the mean- 
square end-to-end distance between monomers n and m, 
([r„ — r m ] 2 ). For ([r„ — r„J 2 ) the perturbation calculation 
has been carried out. The result for infinitely long chains 
reads (s = n — m) [TB] 




Zs TDT ■ (11) 



We see that v(r) corrects the result for ideal chains, Rl(s) — 
b 2 s, in two ways. It increases the effective bond length 
from b to 6 C , as already predicted by Edwards [2J, and 
leads to a swelling of the internal distance between two 
monomers due to the second term — c c *Js. (We refer to c G 
as 'swelling factor' in the following.) This swelling is re- 
sponsable for the deviations of C pp from the ideal behavior 
(cf. Fig. 2]), as we will demonstrate now. 

i?c(s) depends only on the absolute value of the curvi- 
linear distance s. Thus, we first rewrite Eq. Q as 



C, 



1 f N , f" , a/, npir 
/ dn dm R e {\n-m\) cos — cos 

«/ " 

1 + (-1)P+1 



r-JV 



2N 2 



N 



dn / dm R c (n — m) x 
Jo 

npir mqir 

cos cos 

N N 



mqn 
N 



(12) 



and then insert Eq. (fTT]) into Eq. (fl~2"]) . The first term of 
Eq. pT]) will yield the result for ideal chains, i.e. Eq. © 
in the limit p/N <C 1. The second term of Eq. (fTT]) will 
modify Eq. |5]) in two respects: (i) it provides a correction 
to ideal chain behavior (i.e. to the diagonal terms of C pq ), 
(ii) the Rouse modes are not diagonal any longer. 

To substantiate these expectations we have to calculate 
the integral: 



In 



N 

dn 
Jo 



npn mqn 

dm\/n — mcos cos . (Id) 

N N y ' 



Substituting s = n — m and integration by parts gives 
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where S(y/x) and C(y/x) are the Fresnel integrals 



/2tt 
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/2tt 



smy 



cosy 



(15) 



For the term in curly braces in Eq. (|14p we may use the 
expression 



sqn 
cos — — b 

N 



sqn 
~N~ 



• sq-K 
sin C 

N 



sqn sqn 

C08 —- sm w 



sqn\ 



dt c V T * cos i 2 



(16) 



This expression is helpful because the underlined term 
vanishes upon integration over s, whereas the other two 
terms provide the expected corrections to ideal behavior: 
the cosine term amends the diagonal components and the 
third term of Eq. (fT6|) makes the Rouse modes nondiago- 
nal. 

Putting these results together we find (for p, q > 0) 



a 



i hi 



" 2^ N 
2n 2 N y^8 V q 



(ideal) 



[1 + (-1)p+V8 w 

°P1 -3/2, 



P 



dx cos x 



dt c~ 



(17) 



If we now approximate the integral of the underlined term 
by replacing the upper bound pn by oo, we obtain 



dxcosx / die vp xt cost 2 



q/p 



8 [(?/p) + l][Vg7P + l] 
q~ — yr 



(18) 



Table 1. Survey of BSM and BFM parameters for systems 
with variable monomer overlap. E is the energy penalty for full 
monomer overlap, g and b a (g) denote respectively the number 
of monomers in the blob and the statistical segment length, ci 
is the empirical swelling factor discussed in [16]. For the BSM 
we take c\ = 1.3c e and for the BFM, ci = c c [c c is defined in 
Eq. JTT]|]. 
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numerical precision and perhaps refinement of the theory. 
Two refinements are envisageable: It is possible to relax 
the approximation l p = oo' in the underlined integral of 
Eq. (|17p and along with that, to also account for finite- iV 
effects (the latter has been done before in the discussion of 
the form factor [18]). We plan to work on both aspects — 
improvement of numerical precision in the simulation and 
refinement of the theory — in the future. 

Here we continue with the approximation l p = oo' and 
examine its consequences. For the diagonal Rouse modes 
we find 



a 



pp 



Nbl 
2{-Kpf 



1 - 



(for p > 0) . (19) 



and the underlined term vanishes in the limit p — * oo. 

This approximation suggests the Rouse modes to re- 
main essentially diagonal despite of excluded volume inter- 
actions. Indeed, numerical analysis of the theoretical result 
and inspection of the simulation data reveal that the non- 
diagonal terms C pq are smaller than the self-correlation 

C pp by at least two orders of magnitude. The small values 5 Chains with variable monomer overlap 



This result is compared to the simulation data of the 
BSM and the BFM in Fig. [4] For this comparison we 
used the empirical swelling factors c\ determined earlier 
in an analysis of the segmental size distribution in poly- 
mer melts |16j . This analysis revealed that c\ agrees with 
the theoretical value c e for the BFM, while for the BFM, 
a slightly larger value had to be used, c\ = 1.3c c (Ta- 
bic [TJ. Here we adopt these values so that the compari- 
son between theory and simulation contains no further ad- 
justable parameter. Figure[2]demonstrates that the theory 
can indeed account for the systematic depression of C pp 
below the ideal asymptote, if p/N < 0.1. On the other 
hand, the behavior of C pp for larger values of p/N and 
in particular the upturn for p — > cannot be described 
because such large modes should be strongly influenced 
by the microstructure of the polymer model, which is not 
treated correctly by the present theory. 



of C pq make it difficult to separate signal from noise in the 
simulation. However, the results are still indicative of the 
parity property predicted by Eq. (fT7|) : modes with 'p+q = 
odd' vanish, while those with l p+q = even' are finite. This 
parity is a consequence of the translational invariance — 
R e depends only on \s\ — of the chain which we assumed 
to be infinitely long in Eq. ([11"]) . A more quantitative anal- 
ysis of the cross correlations requires improvement of the 



In section [3] we found that phantom chain models are not 
capable of describing conformational properties of flexible 
polymers in a melt by comparing simulation data for ex- 
cluded volume chains to phantom chain calculations. The 
simulation provides a further means to support this con- 
clusion. One can vary the strength of the repulsive inter- 
action between two monomers. 
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Fig. 5. Cpp versus p/N for hard (E = oo) and soft (finite 
values for E) monomers from BSM simulations for N = 256. 
The smaller the value of the overlap energy E, the softer the 
monomers (see Sec. [2}. 



Figure [5] exemplifies the results of such an analysis for 
the BSM. The figure depicts C pp for chains with full ex- 
cluded volume (E — oo) and with finite energy penalty 
E for monomer overlap. Obviously, the p-dependence of 
Cpp qualitatively changes with E. As the monomers be- 
come softer with decreasing E, the curvature of C pp for 
P/N < 0.1, observed for the excluded volume chains, pro- 
gressively vanishes. This clearly demonstrates that this 
feature is determined by the strength of monomer repul- 
sion in the melt. It is possible to rationalize these findings 
by an extension of the perturbation theory of Section [5] to 
soft monomers. 



5.1 Internal distances for finite monomer overlap 

To calculate C pp we begin, as before, by determining R 2 (s) 
first. Quite generally, R 2 {s) i s related to the bond corre- 
lation function 



P\{\n — m\) 



dr(n) dr(m) 



dn 



dm 



(20) 



by 



dn / dm 
Jo 



dr(n) dr(m) 



dn 



dm 



2l 2 







[ dt 


1 - - 


la 


s. 



Pl(t) 



(21) 



where I is the average bond length. For monomers with 
different degree of softness the expression for Pi(s) has 
recently been obtained from first-order perturbation the- 
ory. It reads [21] 



Prix) 
where 



4V7TC P (ff) 
fl 3/2 



1 



- V2e 2x erfc(V2^) 



(22) 
(23) 



Ce(g) is given by Eq. (fTTj) with b c = b c (g), 'erfc' is the 
complementary error function, x = s/g, and g is the num- 
ber of monomers in a blob. (We will comment on the blob 
after Eq. (|2l)|) .) g is related to the correlation length £ and 
the effective bond length b of an unperturbed ideal chain 
via f21] 

r I = — = S(q 0) , (24) 



1 

vp 



where S (q) is the collective structure factor of the polymer 
melt. Inserting Eq. (|22|) into Eq. (j2l"1) we obtain 



R 2 c (s,g) 
b 2 c (g)s 



\[Q 



1 



7r 1 r 
8 xl 



1 - e 2x erfc (V2x 



(25) 



Equation (|25| has the following asymptotic behavior 

R 2 c (s,g)' 



c c (g) 



1 



bl{g)i 



i 



2 5/2 



3tt 



172- 



for x -C 1 , 
for x > 1 . 



(26) 



Here we return to the comment mentioned above. Soft- 
ening of the monomer repulsion introduces the blob size g 
as a new parameter, and the situation becomes similar to 
that of semidilute solutions in good solvent. Inside the blob 
chain segments behave as if they were in dilute solution. 
They do not see their neighbors and are slightly swollen. 
Since the monomer-monomer repulsion is weak, we find 
for iCl that the swelling takes the form of a first-order 
Fixman expansion, familiar from the study of excluded 
volume effects in dilute solutions close to the theta point 
2.3 . For large x — that is, for chain segments much bigger 
than the blob size — Eq. (|25|) gives back Eq. (jTTJ) . On such 
large scales the polymer system behaves like a dense melt 
of blob chains which repel each other. Hence, the swelling 
of internal distances is the same as that of chains whose 
monomers have full excluded interaction. 

To compare simulation and theory we have to deter- 
mine the swelling factor and the blob size. The blob size 
was obtained from the low-g limit of S(q) [cf. Eq. (|24[) ] . For 
large E this limit can be read off reliably from the simula- 
tion data. For the weakest energy penalties, however, S(q) 
does not reach a plateau for the smallest q- values studied. 
Here we determined g by fitting S(q) to the RPA formula 
Siqy 1 = g^ 1 + b 2 c (g)q 2 /12 0. The resulting blob sizes 
have thus larger error bars than those for large E. The fit 
to the RPA formula also yields an estimate for b c (g). This 
estimate can be crosschecked and optimized when deter- 
mining the swelling factors. For the swelling factors we 
adopt the result found previously — that is, c\ = c e for the 
BFM and c\ = 1.3c e for the BSM — but allow c\ to depend 
on g via b e . For small g, i.e. weak to vanishing monomer 
overlap, b 2 {g) may be obtained fairly reliably by fitting 
the asymptotic behavior for x 1 to the simulation data 
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=s/g 



Fig. 6. Scaling plot of Rc(s), as suggested by Eq. (|25|) . for 
monomers of different softness. The upper panel depicts BSM 
data for N = 256, the lower panel shows BFM data for N = 
2048. The dashed line indicates Eq. (f25|) . The solid lines show 
the asymptotic behavior from Eq. (|26[) for small and large x = 
s/g. The values used to scale the axes may be found in Table [T] 




x=gp/N 

Fig. 7. Scaling plot of C pp , as suggested by Eq. (|27|l . for 
monomers of different softness. The upper panel depicts BSM 
data for N = 256, the lower panel shows BFM data for 
N = 2048. The dashed line indicates Eq. p?). The solid lines 
show the asymptotic behavior from Eq. ([28} for small and large 
x = gp/N. The values to scale the axes are given in Table [1] 



for i?c(s, g) [16]. For large g, this fitting procedure is more 
problematic because one has to choose the fit interval, and 
it is hard to find an extended regime (of intermediate x val- 
ues) where the simulation and theoretical curves have the 
same shape. The results for b^(g) were thus obtained by a 
two-step procedure: First, Eq. ([25]) was fitted to Rl(s,g). 
Then, the fit result for b%(g) was optimized in such a way 
that it yields a good data collapse in Fig. [6] and that the 
resulting master curve is close to the theoretical predic- 
tion, Eq. The so-obtained values for g and b^(g) are 
collected in Table [TJ 

In Fig. [6] we compare Eqs. (|25|) and (|26|) with the nu- 
merical results from the BSM and the BFM for various 
values of the energy penalty E. The figure demonstrates 
that there is good agreement between theory and simula- 
tion, perhaps with the exception of the smallest E- values. 
However, here the numerical uncertainties for g and b e are 
largest. 



5.2 Rouse modes for finite monomer overlap 



We find for the diagonal elements of the Rouse mode ma- 
trix the following expression (again in the limit p — > oo, 
see approximation in Eq. (|18[) for comparison) 



2(?rp) 2 r _ , _ /tt c e (g) 



1 



1 



2 ^5 ll + 4/(7rx) 2 

, 7T 1 1 

1 - v 7 ™ + 2 x 



(27) 



where x = gp/N . This equation has the following asymp- 
totic behavior 



c c (g) 



Mr _1 



V8 



- 1 + l/y^rir 



for x -C 1 , 



for x 3> 1 . 



(28) 



The result for the mean-square internal end-to-end dis- 
tance may be inserted into Eq. (J7|) to derive a scaling pre- 
diction for the correlation function of the Rouse modes. 



Equations (J27J) and (gSJ are compared to BSM and BFM 
data in Fig. [7J The agreement between theory and simu- 
lation is of similar quality as in Fig. [6] 
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6 Summary 

The simulation models studied in this work have very flex- 
ible chains in common. As chain stiffness effects are re- 
duced to a minimum — they only come in due to the avoid- 
ance of immediate chain backfolding — one should expect 
ideal chain behavior to appear clearly. 

We demonstrated that, even under these favorable cir- 
cumstances, the assumption of chain ideality on all length 
scales down to the monomer size does not hold. The cor- 
relation function of the Rouse modes displays systematic 
deviations from the (p/N)~ 2 scaling expected for ideal 
chains. Our analysis suggests that these deviations may 
be traced back to the fact that repulsive interactions be- 
tween chain segments in the melt are not fully screened. 
For chain segments of size s 3> 1 there is an entropic 
penalty ~ l/v 7 * f° r bringing two segments together [IB] . 
This penalty swells the segments and causes systematic 
deviations from Flory's ideality hypothesis. The picture of 
independent chains with random-walk-like conformations 
is thus not acceptable for polymer melts. Multi-chain ef- 
fects reflecting the interplay of chain connectivity and melt 
incompressibility should be taken into account. 

In the present work we discussed the impact of these 
multi-chain effects on static chain properties. It is, how- 
ever, natural to expect that they will also affect the poly- 
mer dynamics. An analysis of this influence is underway. 

We are indebted to S. Obukhov for valuable discussions and to 
the IDRIS (Orsay) for a generous grant of computer time. Fi- 
nancial support by the IUF, the ESF STIPOMAT programme, 
and the DFG (grant number KR 2854/1-1) is gratefully ac- 
knowledged. 
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